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Abstract 

We extend the findings and analyses of our two recent studies (Phys. Rev. A 75, 032326 



[2007] and arXiv:0704.3723) by, first, obtaining numerical estimates of the separability function 
based on the (Euclidean, flat) Hilbert-Schmidt (HS) metric for the 27-dimensional convex set of 
quaternionic two-qubit systems. The estimated function appears to be strongly consistent with our 
previously-formulated Dyson-index {(5 = 1,2,4) ansatz, dictating that the quaternionic {(3 = 4) 
separability function should be exactly proportional to the square of the separability function for 
the 15-dimensional convex set of two-qubit complex {(5 = 2) systems, as well as the fourth power 
of the separability function for the 9-dimensional convex set of two-qubit real (f3 = 1) systems. 
In particular, we conclude that S^ t (fi) = (|[) 2 ((3 - u 2 )u^j = (S^ plex {u)) 2 , < u < 1. Here, 
= \ q 11 ^ 44 ) where p is a 4 x 4 two-qubit density matrix. We can, thus, supplement (and fortify) 
our previous assertion that the HS separability probability of the two-qubit complex states is 
w 0.242424, by claiming that its quaternionic counterpart is 9 7 3 2 6 4 2 4 3 2 9 9 7 4 2 4 5 w 0.0773765. We also 
comment on and analyze the odd (3 = 1 and 3 cases. 
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For several years now, elaborating upon an idea proposed in we have been pursuing 



the problem of deriving (hypothetically exact) formulas for the proportion of states o 
}it and qubit-qutrit systems that are separable (classically-correlated) in nature 
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9J. Of course, any such proportions will critically depend upon the measure 



that is placed upon the quantum systems. In particular, we have — in analogy to Bayesian 
analyses, in which the volume element of the Fisher information metric for a parameterized 



family of probability distributions is utilized as a measure ( "Jeffreys' prior" ) lOj — principally 
employed the volume elements of the well-studied (Euclidean, flat) Hilbert-Schmidt (HS) and 
Bures (minimal monotone) metrics (as well as a number of other [non- minimal] monotone 
metrics 6j). 

Zyczkowski and Sommers 



11 



121 ] have, using methods of random matrix theory [13( (in 
particular, the Laguerre ensemble), obtained formulas, general for all n, for the HS and 
Bures total volumes (and hyperareas) ofnxn (real and complex) quantum systems. Up 
to normalization factors, the HS total volume formulas were also found by Andai [141 ]. in 
a rather different analytical framework, using a number of (spherical and beta) integral 
identities and positivity (Sylvester) conditions. (He also obtained formulas — general for any 
monotone metric [including the Bures] — for the volume of one-qubit [n = 2] states [l4j, sec. 
4]-) 

Additionally, Andai did specifically study the HS quaternionic case. He derived the HS 
total volume for n x n quaternionic systems p. 13646], 

v hs _ (2n-2)!7r n ~ n n n-2/ 2 - v m 
giving us for the two-qubit (n = 4) case of specific interest here, the 2 7- dimensional volume, 
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2.93352 ■ 10~ 18 . (2) 



(In the analytical setting employed by Zyczkowski and Sommers llj, this volume would 



appear as 2 13 times as large [14J, p. 13647].) If one then possessed a companion volume for- 
mula for the separable subset, one could immediately compute the HS two-qubit quaternionic 
separability probability by taking the ratio of the two volumes. 

One analytical approach to the separable volume/probabrlity option that has recently 
proved to be productive [15j — particularly, in the case of the Hilbert-Schmidt (HS) metric (cf. 



16j) — makes fundamental use of a form of density matrix parameterization first proposed 



by Bloore 17J. (This methodology can be seen to be strongly related to the very common 
and long-standing use of correlation matrices in statistics and its many fields of application 
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20|.) 



In the Bloore parameterization, one simply represents an off-diagonal zfentr y of a density 

23j in 



matrix p, as py 



V PaPjj w ij ; 



where Wij might be real, complex or quaternionic [2JJ 
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nature. The particular attraction of the Bloore scheme, in terms of the separability problem 



in which we are interested, is that one can (in the two-qubit case^ 
known Peres-Horodecki separability (positive-partial-transpose) test 
ratio p = 



implement the well- 



24 . |25| using only the 



P11P44 , rather than the four (three independent) diagonal entries of p individually 
B, eq. (7)] V eq. (5)]. 

Utilizing the Bloore parameterization, we have, accordingly, been able to reduce the prob- 
lem of computing the desired HS volumes of two-qubit separable states to the computations 
of one-dimensional integrals over p e [0,oo]. The associated integrands are the products 
of two functions, one a readily determined jacobian function J(p) (corresponding, first, to 
the transformation to the Bloore variables Wij and, then, to p) and the other, the more 
problematical (what we have termed) separability function S HS (p) 0, eqs. (8), (9)]. (In the 
qubit- qutrit case, two ratios, p\ and p2, are required to express the separability conditions, 
but analytically the corresponding separability functions also appear to be univariate in 
nature, being simply functions of p\ or p 2 , or the product p\p 2 [la, sec. III].) 

In our extensive numerical (quasi-Monte Carlo integration) investigation 9] of the 9- 
dimensional and 15-dimensional convex sets of real and complex 4x4 density matri- 
ces, we had formulated ansatze for the two associated separability functions (S^(p) and 
<S^ pleT (fj, )\ p roposing that they were proportional to certain (independent) incomplete beta 
functions 261 ]. 
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for particular values of a and b. However, in the subsequent study [la], we were led to 
somewhat modify these ansatze, in light of multitudinous exact /ower-dimensional results. 
Since these further results clearly manifested patterns fully consistent with the Dyson index 
("repulsion exponent") pattern (/3 = 1, 2, 4) of random matrix theory |27|, we proposed that, 
in the (full 9-dimensional) real case, the separability function was proportional to a specific 



3 



incomplete beta function (a = |, b = 2), 

SZZWBpb?, ^,2)^(3-/i 2 )/i (4) 

and in the complex case, proportional, not to an independent function, but simply to the 
square of (These proposals are strongly consistent 15|, Fig. 4] with the numerical 

results generated in Qj.) This chain of reasoning, then, immediately compels one to the 
further proposition that the separability function in the quaternionic case is exactly propor- 
tional to the fourth power of that for the real case (and, obviously, the square of that for 
the complex case). It is that specific proposition we will, first, seek to evaluate here. 

We, thus, hope thereby to further test the validity of our Dyson-index ansatz, first ad- 
vanced in [15J, as well as possibly develop an enlarged perspective on the still not yet fully 
resolved problem of the HS separability probabilities in all three (real, complex and quater- 
nionic) cases. (In , we proposed, combining numerical and theoretical arguments, that in 
the real two-qubit case, the HS separability probability is ^, and in the complex two-qubit 
case, Jj. The arguments, thusly, employed in 
formal demonstration.) 

Due to the "curse of dimensionality" 



151 ]. however, do not yet rise to the level of a 
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291 ]. we must anticipate that for the same 



number of sample ("low-discrepancy" Tezuka-Faure 30|, |3l||) points generated in the quasi- 
Monte Carlo integration procedure employed in 9J and here, our numerical estimates of 
the quaternionic separability function will be less precise than the estimates were for the 
complex, and a fortiori, real cases. (An interesting, sophisticated alternative approach 



32] 



to computing the volume of convex bodies involves a variant of simulated annealing 
(cf. [33!]), and allows one — unlike the Tezuka-Faure approach, we have so far employed — to 
establish confidence intervals for estimates.) 

Our first extensive numerical analysis here involved the generation of sixty-four million 24- 
dimensional Tezuka-Faure points, all situated in the 24-dimensional unit hypercube [0, 1] . 
(The three independent diagonal entries of the density matrix p — being incorporated into 
the jacobian 3{\i) — are irrelevant at this stage of the calculations of S^ at {p). The 24 [off- 
diagonal] Bloore variables had been transformed so that each ranged over the unit interval 
[0,1]. The computations were done over several weeks, using compiled Mathematica code, 
on a MacMini workstation.) 

Of the sixty-four million sample points generated, 7,583,161, approximately 12%, corre- 
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sponded to possible 4x4 quaternionic density matrices — satisfying nonnegativity require- 
ments. For each of these feasible points, we evaluated whether or not the Peres- Horodecki 
positive-partial-transpose separability test was satisfied for 2,001 equally-spaced values of 

[0,1]. 

Here, we encounter another computational "curse" , in addition to that already mentioned 
pertaining to the high-dimensionality of our problem, and also the infeasibility of most (88%) 



of the sampled Tezuka-Faure points. In the standard manner 
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eq. (5.1.4)] Q, p. 495] 
34|, eq. (17)] 35|, sec. II], making use of the Pauli matrices, we transform the 4x4 
quaternionic density matrices — and their partial transposes — into 8x8 density matrices 
with [only] complex entries. Therefore, given a feasible 24-dimensional point, we have to 
check for each of the 2,001 values of ft, an 8 x 8 matrix for nonnegativity, rather than a 
4x4 one, as was done in both the real and complex two-qubit cases. In all three of these 
cases, we found that it would be incorrect to simply assume — which would, of course, speed 
computations — that if the separability test is passed for a certain /z , it will also be passed 
for all ft lying between fto and 1. This phenomenon reflects the intricate (quartic both in ft 
and in the Bloore variables Wi/s, in the real and complex cases) nature of the polynomial 
separability constraints |9j, eq. (7)] [l5j, eq. (5)]. 

In Fig. [1] we show the estimate we, thus, were able to obtain of the two-qubit quaternionic 
separability function <Sg U a t {fi), in its normalized form. (Around fi — 1, one must have the 
evident symmetrical relation S HS (fi) = S HS (j^).) Accompanying our estimate in the plot 
is the (well-fitting) hypothetical true form (according with our Dyson- index ansatz [l5|]) of 
the HS two-qubit separability function, that is, the fourth power, ^|(3 — fi 2 )fi) , of the 
normalized form of <S^(/i). 

For the specific, important value of fi — 1, the ratio (Ri) of the 24-dimensional HS 
measure (m sep = R™ umer ^ assigned in our estimation procedure to separable density matrices 
to the total 24-dimensional HS measure (m tot = Rdenom^ a u ^ed to all (separable and 
nonseparable) density matrices is Ri = 0.123328. The exact value of m sep is, of course, to 
begin here, unknown, being a principal desideratum of our investigation. On the other hand, 
we can directly deduce that m to t = Rf enom = 77 jH 00 ~ 0.118862 — our sample estimate being 
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FIG. 1: Estimate — based on 64,000,000 sampled 24-dimensional points — of the normalized form 
of the two-qubit quaternionic separability function, along with its (well-fitting) hypothetical true 
form, the fourth power of the normalized form of S^^fi), that is, ^|(3 — u 2 )u^j 

0.115845 — by dividing the two-qubit HS quaternionic 2 7- dimensional volume (j2j) by 

r 1 rf^ + i) 4 1 

R denom = 2 / j ' (n)du = — ^ k~ = « 2.46801 • 10" 17 , = 4. 

2 ^ uquat^j f" y(6[3 + 4) 40518448303132800 

(5) 

Here, J qua t(p) i s the quaternionic jacobian function (Fig. [2]), obtained by transforming the 

3/3 

quaternionic Bloore jacobian \ p\\p-nPm*(y — Pu — P22 — P33) j 2 , /3 = 4, to the p variable by 
replacing, say P33 byu, and integrating out p\\ and pii- (We had presented plots of J re ai(p) 
and Jcompiex{.p) m [91, Figs. 1, 2], and observed apparently highly oscillatory behavior in 
both functions in the vicinity of p = 1. However, a referee of JJjJ informed us that this was 
simply an artifact of using standard machine precision, and that with sufficiently enhanced 
precision, the oscillations could be seen to be, in fact, illusory.) We can obtain an estimate 



of the two-qubit quaternionic separability probability P!r e p/ quat by multiplying the ratio R\ 
by a second ratio i?2- The denominator of R2 has already been given (jSJ). The numerator 
of i?2 is the specific value 



r 1 f \ n4 5Q8Q 

R numer = 2 / J Jn) ( - (3 - n 2 ) n) dp = « 1.67128 -10" 

2 y Q u q uatw)y 2 \ p jh-j p 358347086242825680000 



17 



(6) 



where, to obtain the integrand, we have multiplied (in line with our basic [Bloore- 
parameterization] approach to the separability probability question) the quaternionic jaco- 
bian function by the (normalized) putative form of the two-qubit quaternionic separability 
function. (Note the use of the (3 = 4 exponent.) 
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FIG. 2: The univariate quaternionic jacobian function J qU at( a ) 



The counterpart of R™ umer in the 9-dimensional real case is 



i 



151200 



and in the 15- 



dimensional complex case 
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' 99891792000 



We now note that 



99891792000 




£(16) 

r(7) 



(7) 



is the coefficient of /i 2 in ll!L^ 1 (/x) and 



151200 



75600 plays the exact same role in 6!L 4 (/i), 



where is a generalized (a = 4) Laguerre polynomial (see sequences A062260 and 

A062140 in the The On-Line Encylopaedia of Integer Sequences). (Also, as regards the 
denominator of ® , 35834 g"^^^ 80000 = 99 89 1 79 2 00.) Zyczkowski and Sommers had made 

e in deriving the HS and Bures volumes and hyperareas of 



use of the Laguerre ensemb 
n-level quantum systems 11 



121 ]. Generalized (associated/Sonine) Laguerre polynomials 



"Laguerre functions") have been employed in, in another important quantum-information 



context, in proofs of Page's conjecture on the average entropy of a subsystem 
We, thus, have, for our two-qubit quaternionic case, that 



37j.) 



R-2 



R numer 125769 



TDdenom 



185725 



0.677179. 



(8) 
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(The real counterpart of R2 is ~ 0.76854, and the complex one, || ~ 0.717172. Ad- 
ditionally, we computed that the corresponding "truncated" quaternionic 38( ratio — when 



one of the four quaternionic parameters is set to zero, that is the Dyson-index case (3 = 3— 
* s 106376244975^ ~ 0.692379. Thus, we see that these four important ratios monotonically 
decrease as j3 increases, and also, significantly, that the two ratios for odd values of f3 differ 
qualitatively — both having 7r 2 in their denominators — from those two for even (3.) 

Our quasi-Monte Carlo (preliminary) estimate of the two-qubit quaternionic separability 
probability is, then, 

Ps H ep/ q uat « #1#2 = 0.0813594. (9) 

Multiplying the total volume of the 27-dimensional convex set of two-qubit quaternionic 
states, given in the framework of Andai [yj] by ([2]), by this result ([9]), we obtain the two- 
qubit quaternionic separable volume estimate V^ quat ~ 2.38775 ■ 10~ 19 . 

Our 24-dimensional quasi-Monte Carlo integration procedure leads to a derived estimate 
of (the total 27-dimensional volume) V^ t , that was somewhat smaller, 2.85906 • 10~ 18 , 
than the 2.93352 • 10~ 18 given by (T5]). Although rather satisfying, this was sufficiently 
imprecise to discourage us from attempting to "guestimate" the (all-important) constant 
(i?i) by which to multiply the putative normalized form, (|(3 — /i 2 )/i) 4 , of the quaternionic 
separability function in ([6]) in order to yield the true separable volume. In our previous 



study 



15l . sec. IX. A], we presented certain plausibility arguments to the effect that the 



corresponding constant in the 9-dimensional real case might be = (^f-)/( 51 2 2 ^ ), and 
if = ( 2 639 6 )/(^ff~) i n the 15-dimensional complex case. (This leads — multiplying by the 
corresponding iVs, and || — to separability probabilities of ^ and Jjj, respectively.) 

In light of such imprecision, we undertook a supplementary analysis, in which, instead 
of examining each feasible 24-dimensional point for 2,001 possible values of fi, with respect 
to separability or not, we simply used /i = 1. This, of course, allows us to significantly 
increase the number of 24-dimensional Tezuka-Faure points generated from the 64,000,000 
so far employed. 

We, thusly, generated 1,360,000,000 points, finding that we obtained a remarkably good 
fit to the important ratio Ri of the 24-dimensional measure, at \i — 1, assigned to the 

12 

separable two-qubit quaternionic density matrices to the measure (known to be 77 y 60Q0 ) by 
setting Ri = (ff) 2 ~ 0.114263 (our sample estimate of this quantity being 0.114262). This 
is exactly the square of the corresponding ratio lj we had conjectured (based on extensive 

8 



numerical and theoretical evidence) for the full (15-dimensional) complex two-qubit case in 

Q. 

Under this hypothesis on Ri, we have the ensuing string of relationships 

Oa0 = (^) 2 (^( 3 -^) 4 = Q 2 (( 3 -^) 4= { s coi Ple M)\ (io) 



with (as already advanced in 151 ]) 

94/1 \ 2 fi / \2 

S^ P ieM = ji(^-n 2 )n) =^((3-^) . (11) 

Then, using our knowledge of the complementary ratio R2, given in (jSJ), we obtain 

p 5Luat = R1R2 = 72442944 ~ 0.0773765, (12) 
sep/guat 936239725 V 1 



as well as — in the framework of Andai 



ij— that 



KfiLw = « 2.26986 • 10~ 19 . (13) 

sep/gnat 24386773433626137413880000000 v ; 

For possible further insight into the HS two-qubit separability probability question, we 
undertook a parallel quasi- Monte Carlo (Tezuka-Faure) integration (setting \i — 1) for the 
truncated quaternionic case (/3 = 3), in which one of the four quaternionic parameters is 
set to zero. Although there was no corresponding formula for the HS total volume for this 



scenario given m 



141 ]. upon request, A. Andai kindly derived the result 
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Vt ™ nc = 384458588946432000 ~ 2,43584 ' (14) 
(In fact, Andai was able to derive one simple overall comprehensive formula — which we leave 
for him to publish — yielding the total HS volumes for all n x n systems and Dyson indices 
(3.) Let us, further, note that Andai obtains the result (CHI) as the product of three factors, 
v trunc = tti^tts, where 

128tt 8 128 189tt 2 

7Ti = : 7T9 = ; 7T3 = . 15 

105 893025 12696335643836416 v ; 

Now, we will simply assume — in line with our basic Dyson- index ansatz, substantially 
supported in [3| and above — that the corresponding separability function is of the form 

^ 5 nc (/i)cx((3-/i 2 )/if, (3 = 3. (16) 

(Of course, one should ideally test this specific application of the ansatz too, perhaps in the 
manner we have examined the (3 = 4 instance above [Fig. [[].) 



We were somewhat perplexed, however, by the results of our quasi-Monte Carlo integra- 
tion procedure, conducted in the 18-dimensional space of off-diagonal entries of the truncated 
quaterionic density matrix p. Though, we anticipated (from our previous extensive numeri- 
cal experience here and elsewhere) that our estimate of the associated 18-dimensional volume 
would be, at least, within a few percentage points of 7ri7r 2 = 93767625 ~ 1-65793, the estimate 
was, in fact, close to 0.967 (1, thus, falling within the possible margin of error). Assuming 
the correctness of the analysis of Andai, which we have no other reason to doubt, the only 
possible explanations seemed to be that we had committed some programming error (which 
we were unable to discern) or that we had some conceptual misunderstanding regarding the 
analysis of truncated quaternions. (Let us note that we do convert the 4x4 density matrix 
to 8 x 8 [complex] form 22j, p. 495] 34j, eq. (17)] 35|, sec. II], while it appears that Andai 
does not directly employ such a transformation in his derivations.) 

In any case, we did devote considerable computing time to the (3 = 3 problem (generating 
1,180,000,000 18-dimensional Tezuka-Faure points), with the hope being that if we were in 
some way in error, the error would be an unbiased one, and that the all-important ratio of 
volumes would be unaffected. 

Proceeding thusly, our best estimate (not making use of the Andai result ffl4|) for the 
present) of the HS separability probability was 0.193006. One interesting possible candidate 
exact value is, then, ||| = 3 ,^, 17 ~ 0.193062. (Note the presence of 128 in the numerators, 
also, of both factors 7i"i and 7^.) This would give us a counterpart [(3 = 3] value for the ratio 



R 2 of 1 !S?6i6 ~ °- 278838 - In U3|, we had asserted that, in the other odd (3 = 1 case, the 
counterpart of R2 was 1 2 3 1 5 ^ 2 ~ 0.612315. (Multiplying this by gave us the conjectured 
HS real two-qub it separability probability of ^.) 

So, let us say in conclusion, that although we believe we have successfully resolved — 
though still far from having formal proofs — the two-qubit Hilbert-Schmidt separability prob- 
ability question for the (3 = 2 and 4 (complex and quaternionic) cases, the odd ((3 = 1, 3) 
cases, in particular (3 = 3, appear still to be somewhat more problematical. 
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